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The hopping motion of lattice gases through potentials 
without mirror-reflection symmetry is investigated under var- 
ious bias conditions. The model of 2 particles on a ring with 
4 sites is solved explicitly; the resulting current in a sawtooth 
potential is discussed. The current of lattice gases in extended 
systems consisting of periodic repetitions of segments with 
sawtooth potentials is studied for different concentrations and 
values of the bias. Rectification effects are observed, similar to 
the single-particle case. A mean-field approximation for the 
current in the case of strong bias acting against the highest 
barriers in the system is made and compared with numerical 
simulations. The particle-vacancy symmetry of the model is 
discussed. 



I. INTRODUCTION 

The motion of particles in potentials that do not have 
mirror reflection symmetry has attracted much attention 
in the last years for several reasons. The interest extends 
from fundamental problems concerning the validity of the 
Second Law of Thermodynamics |^J^ to applications in 
biological and chemical systems |8j, as well as for 

solid-state devices |9|[Tot. Major efforts have been de- 
voted to an understanding of molecular motors, where 
proteins move in nonsymmetric potentials under the in- 
fluence of stochastic and/or other forces. One specific 
observation for transport in nonsymmetric potentials is 
the possibility of rectification effects if the forces on the 
particles are beyond the regime where linear-response 
theory is applicable Rectification effects have been 
discussed in continuous as well as in hopping sys- 
tems 1^,0 . If applications of effects of particle motion in 
nonsymmetric potentials are envisaged then the question 
arises as to the influence of many-particle effects. The 
limit of single-particle motion is rarely realized; in real 
systems many particles are present that compete about 
the sites that can be occupied. Many-particle effects have 
been studied in continuous nonsymmetric periocic poten- 
tials in Ref. ||l^, where interesting dependencies of the 
current on particle concentration and size were found. In 
this paper we will investigate hopping motion of lattice- 
gas particles in nonsymmetric hopping potentials under 
the influence of strong bias. We utilize the simple site 
exclusion model where multiple occupancy of sites is ex- 
cluded and direct our attention to nonlinear effects on 



the particle current. 

The stationary current of a single particle performing 
a hopping motion in a nonsymmetric potential under an 
arbitrary bias is known exactly The calculation of 
the stationary current of site-exclusion lattice gases in 
nonsymmetric potentials that lead to rectification effects 
in the single-particle case is a difficult problem. Exten- 
sive work has been devoted to the asymmetric site ex- 
clusion process including the totally asymmetric site ex- 
clusion process (TASEP) where the particles can only 
hop in one direction, corresponding to very strong bias. 
The case of uniform hopping potentials is now well un- 
derstood jl^, but the case of nonuniform potentials is 
not generally solved. Recent work has been devoted to 
the TASEP with disordered potentials For the 

general asymmetric case one has to resort to numerical 
simulations; we are going to present simulation results for 
the stationary current of lattice-gas particles in nonsym- 
metric potentials, for various concentrations and values 
of the bias. 

Nonetheless, some analytical treatment can be given. 
First, the case of very small periodic systems can be 
treated explicitly: the motion of 2 particles on a ring of 
period 4 can be solved by elementary means. Although 
this is a very simple system, conclusions can be drawn in 
the limit of very strong bias that are of interest for the 
totally asymmetric site exclusion process. The nonlinear 
current of site-exclusion lattice gases in extended systems 
with periodic repetitions of nonsymmetric segments can 
be derived in a mean-field approximation for strong bias 
conditions. Interesting symmetry properties have been 
pointed out for the TASEP in disordered hopping poten- 
tials [p^|-pO|. While a particle- vacancy symmetry is also 
present in our model, the case of inversion of the bias 
direction is different here. 

In the following section the hopping motion of 2 par- 
ticles on a ring of period 4 is solved and analyzed. In 
Sec. Ill a mean-field approximation for the stationary 
current of lattice gases under strong bias in nonsymmet- 
ric hopping potentials is presented and compared with 
simulation results in a sawtooth potential. The symme- 
try properties of the model are discussed in Sec. IV and 
concluding remarks are given in Sec. V. 



II. TWO PARTICLES ON RING WITH FOUR 
SITES 
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A. Solution of the stationary master equations 

A very simple yet nontrivial model is given by a ring 
with 4 sites and 2 particles, cf. Fig. |l|. The basic quanti- 
ties for the description of the system are the joint prob- 
abilities P{i,j;t) {i / j) of finding one particle at site i 
and the other particle at site j, at time t, for specified 
initial conditions. Since the particles are considered as 
indistinguishable, P{i,j;t) — P{j,i;t). There are 6 dif- 
ferent joint two-particle probabilities on the ring with 4 
sites (generally L{L — l)/2 on rings with L sites). Higher- 
order joint probabilities do not occur for 2 particles. 

The probabilities P{i\t) of finding a particle at site i 
at time t are given by 



joint probabilities for adjacent sites, e.g. P12, can all be 
expressed by the stationary joint probabilities P13 and 
P24- For instance, the first line of Eq.M) yields 



P12 



1 



72 + ^1 



((53^13 + 74^24) 



(6) 



Three analogous relations follow from (^); they can be 
obtained by cyclically increasing the indices in Eq. (||) . If 
the joint probabilities for adjacent sites are eliminated 
from the stationary master equations, two homogeneous 
equations remain which are equivalent. We write this 
equation as 



p{i;t)^Y.P^^-'r^^) 

For two particles they are normalized to 
4 



(1) 



ailPl3 + 012^24 = 

with the coefficients 

6i6^ 7173 6163 



(7) 



ail 



ai2 



72 + 5i 



^ ^ 7173 ^ 

73 + 52 74 + (^3 71 -f (54 ' 

7274 , 625 a 



This condition implies 

Y.P{i,r,t) = i 



(2) 
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The second equation for P13 and P24 is obtained from 
the normalization condition Eq.(^, after elimination of 
the joint probabilities of adjacent sites. It reads 



The master equations for the joint probabilities are easily 
written down, 
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+ 1. (10) 



Pl3 = 

P24 = 
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an 

aiia22 ~ 012021 



(11) 
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Since the joint probabilities for adjacent sites are ob- 
tained from the P13, P24, and the one-site stationary 
probabilities Pi = P{i;t — > 00) from Eq.(0), Eq. ([lT|) rep- 
resents the complete solution of the stationary problem. 

We derive the stationary current in the system by con- 
sidering the bond connecting sites 1 and 2. The station- 
ary current is given by 



The sum of the 6 master equations leads to the conser- 
vation law 



J = 7l(Pl-Pl2)-<52(P2-Pl2) 



(12) 



-EP(z,j;i)]=0, 



(5) 



i<j 



The joint probabilities in Eq.([T^) ensure exclusion of dou- 
ble occupancy of sites. Using Eq.(|l|) the current is ex- 
pressed in terms of the joint probabilities. 



consistent with the relation (|[) given above. 

We are interested in the stationary solution of the 
system of master equations (^. The stationary values 
P{i,j;t 00) will be denoted by Py. The stationary 



J = 71(^3 + P14) - <52(P23 + P24) 



(13) 



Insertion of the stationary solution for the joint proba- 
bilities gives 
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(71 + 0i){j3 + O2) 

The current may also be derived by considering the other 
bonds of the ring. Two equivalent forms of the current 
result; the second (equivalent) form reads 

(72 + 0l)(74 + O3) 

It can be shown that the current vanishes if the right and 
left transition rates fulfill the following condition, 

71727374 = S16263S4 (16) 
corresponding to a detailed balance relation over the ring. 

B. Solution for the sawtooth potential 

The sawtooth potential including bias on a 4-sitc ring 
is defined by choosing 



71 = 72 = 73 = 74 = bj, 

5i^b-^^\ 

82 = 53 = 84 = 



(17) 



where b represents the bias and 7 < 1 is a constant rep- 
resenting a transition rate to the right in the absence 
of a bias, cf. Fig. 1(b). Note that the right transition 
rates are explicitly multiplied by the bias factor b and 
the left transition rates by 6~^, respectively. Physically, 
b ~ exp{AU/2kBT) where AU represents the potential 
drop between two neighboring sites under the influence 
of the bias. For 6=1 the system satisfies the detailed 
balance condition and the current J vanishes. In what 
follows the current obtained in a system with M particles 
will be denoted as Ja/. 

In Fig. ^ we present a plot of Ji and J2 as functions of 
the bias b for the ring with 4 sites and 7 = exp(— 2). The 
result for the two-particle system was obtained using Eq. 
(|l5|), and for a single-particle system we employed the 
exact solution derived in Ref. We can see that the 
behavior of the currents of one- and two-particle systems 
are qualitatively similar. Of course, the current J2 of 
two particles is larger than the one-particle current Ji. 
The inset shows the behavior of the current for smaller 
bias. The curves for the bias to the right and to the 
left become equal in the limit 6 — > 1, i.e., in the linear- 
response regime, for two particles on the ring with 4 sites, 
and also for one particle on this ring. However, the two- 
particle current is about 17% larger than the one-particle 
current. 

In the case of a strong bias to the right, 6^1, the two- 
particle current J2 differs from Ji by a constant factor. 
For the sawtooth potential this behavior can be under- 
stood as follows. If 6 3> 1, only transitions to the right are 
important, and backward transitions can be neglected. In 



our model the transition rates to the right are all equal, 
7i = 67 for i = 1, ... ,4. Hence for & 3> 1 all station- 
ary site occupation probabilities become equal. Pi = 1/2 
for the 4-site ring and Pi = 2/L for a ring with L sites. 
In the limit of a strong bias to the right all stationary 
joint probabilities also become equal, i.e., V^.j Pij = 1/6 
for the 4-site ring and, generally, Pij = 2/L{L — 1) (see 
Ref. [0). Using expression ( |l2|) we thus expect that for 
&> 1 



J^bj 



L L{L-l) 



(18) 



For L = 4 there is thus J2 = 67/8, which should be com- 
pared with the single-particle current Ji = 67/4. Simi- 
larly, in the general case of an L-site ring we have 



lim 



J2 ^ 2(£-2) 

Ji ^ 



L - 1 



(19) 



For the 4-site ring this limiting behavior can be easily 
derived from the exact formula (|l5|). Actually, for i = 4, 
7 = exp(— 2), and 6 = 10 there is J2/J1 ~ 1.3337, in 
agreement with the above considerations. 

Figure || also shows that J2 becomes almost identical 
to Ji in the case of a strong bias to the left, 6 <C 1. 
To understand this phenomenon assume that 7 1, so 
that 5i <^ 82 = 83 = 84, i.e., site 1 acts as a "bottle- 
neck" . If 6 <C 1 the particles are driven against the high 
barrier at site 1, which has a relatively very small tran- 
sition rate Si to the left. The second particle on site 2 
has to wait until the first particle has jumped over the 
high barrier, and only then can it make an attempt to 
jump over that barrier. Soon after the first particle has 
managed to pass the bottleneck at site 1, the second par- 
ticle will jump from site 2 to 1 and the first particle will 
quickly line up behind the second particle, waiting for it 
to jump over the high barrier. Consequently, the current 
becomes practically equal to that of a single-particle sys- 
tem. It is evident that in the limit of a large bias to the 
left the system behaves as a TASEP on a ring with one 
defect. If the defect is characterized, in a discrete-time 
dynamics, by the transition probability p ^ 1, the cur- 
rent of M particles on a ring with L sites {M < L) will 
approach the one-particle current. The above reasoning 
is confirmed by an explicit calculation of the current J2 
in the limit 6 — s- 0. Using dl3) we conclude, after some 
algebra, that J2 « 2b-^{l + i^){b-i'^ + 2-i~^ + . Since 
for a single-particle system the current Ji, for 6 <C 1, is 
approximately equal 6~^7^(1 -I- 87^^) (see Ref. [|ll|), we 
find that 



lim — = 

b^o Ji 57 



2(l + 7")(l + 37") 
+ 574 + 2 



(20) 



For 7 — > 0, i.e., for a growing asymmetry of the sawtooth 
potential, this limit actually approaches 1. In particu- 
lar, for the value of 7 = exp(— 2) used in Fig. 2 there is 
limbic J2/J1 ~ 1.0005. 
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Note, however, that in contrast to the case 6 ^ 1, for 

6 <C 1 the current depends on the parameter 7 charac- 
terizing the inhomogeneity of the sawtooth potential. In 
particular, for 7=1, which corresponds to a fully ho- 
mogeneous system, \imh-,oJ2/Ji — 4/3. Actually, for 

7 = 1, the ratio J2/J1 equals 4/3 irrespective of the bias 
b (see M). 



III. EXTENDED NONSYMMETRIC 
POTENTIALS 

A. The model 

In this section lattice gases in extended potentials are 
considered that consist of periodic repetitions of nonsym- 
metric segments. First the situation of very strong bias is 
discussed and a mean-field approximation is given for the 
case where the particles experience periodically arranged 
high barriers. The analytical results are then compared 
with numerical simulations of the motion of lattice-gas 
particles in nonsymmetric hopping potentials for different 
concentrations and under various bias conditions. The 
hopping potential that is used in in this section is the 
sawtooth potential as shown in Fig. 1(b), except that 
it is periodically repeated with period L. The nearest- 
neighbor transition rates from site I to I ±1 are Ti i±i. 
As a short notation we use 7; = Tij+i for the "right" and 
61 = Tij-i for the "left" transition rates. Without ad- 
ditional bias, the transition rates between neighbor sites 
fulfill detailed balance. Bias is introduced by multiply- 
ing all right transition rates by 6, 7; — > bji, and all left 
transition rates by b~^, Si — > b~^Si. 

The linear chain on which the model is defined shall 
have N = vL sites where we consider v ^ \ 'm this sec- 
tion. Periodic boundary conditions are introduced and 
the sites are occupied by M particles. The concentration 
is then p = M/N. Multiple occupancy of the sites is ex- 
cluded; no further interactions of the particles are taken 
into account. 



B. Strong bias 

1. The case 6 » 1 

For b ^ I we can apply essentially the same reasoning 
as in the case of the 2-particle system considered in Sec. 
II. In this limit transitions to the left are so rare that 
they can be ignored and the system essentially behaves 
like a TASEP with transitions 7^ = 67, i = 1, . . . , A^. The 
current for such a system reads [ pl| 



J = 67 Af 



N -M 
N - 1 



(21) 



J{p) = bjp{l - p) . (22) 



2. The case 6 < 1 

For 6 <C 1 we can neglect transition rates to the right, 
and so the system behaves like a TASEP with transition 
rates 5i = b"^^'^ if i = 1 (mod L) and 5i = b~^ otherwise. 
If additionally 7 = 1, all (5^ are eq ual to each other and 
the current is given simply by (21). 

A more complicated situation appears for 7 ^ 1, a 
condition which will be assumed henceforth. In this case 
sites i = 1, L + 1, . . . , N — L + 1 act on the flow of particles 
as "bottlenecks", for the mean time necessary to leave 
them is much larger than the time to leave any other 
site. Therefore the system, which consists of v similar 
segments of length L, effectively behaves like a ring made 
up of ly similar "boxes" , each able to contain up to L 
particles. A transition from a segment j to j — 1 occurs 



with a rate 6~ 7 irrespective of the number of particles 
in each of the segments, provided, of course, that there 
is at least one particle in segment j and at most L — 1 
particles in segment j — 1. 

Let Qn denote the probability that in the steady state 
there are n particles in a given segment {n = 0, . . . , L). 
Let Qm.n denote the joint probability of finding, in the 
steady state, < m < L particles at a given segment j 
and < n < L particles at j + 1. Of course Qn and Qm.n 
do not depend on j, and the Qn satisfy 



^Qn = l, 
L 

nQn = Lp. 



(23) 
(24) 



Let us assume a mean-field approximation: Qm.n = 
QmQn- In the stationary state the mean number Mn 
of segments occupied by n particles does not depend on 
time. As the particles hop between segments, Afn can de- 
crease when one of the particles jumps from or to a seg- 
ment occupied by n particles. The corresponding rates 
are Qn(l — Ql) and Qn{^ — Qo), respectively. The num- 
ber of segments containing n particles can also increase 
owing to jumps ending at segments containing n—1 par- 
ticles or originating at segments with n + 1 particles; 
the corresponding transition rates are Q„_i(l — Qo) and 
(5„+i(l — Ql), respectively. Consequently, the appropri- 
ate balance conditions read 



(Q. 



-i)(1-Ql) = il-Qo)iQn-l 

Qi(i-QL)-(i-go)Oo, 

3l(1-0l) = (I-Qo)Ol-i, 



Qn) 



(25) 
(26) 
(27) 



For large system sizes 3> 1 this formula can be rewrit- 
ten as 



where n = 1, . . . , L - 1 in (gS]) and in (|26D and (|23) 
we have taken into account the fact that neither jumps 
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from a segment containing particles nor transitions to 
a segment with L particles are possible. Together with 
( p3| ) and (24) these relations form L + 3 equations for 
L + 1 variables Q„, n — 0, . . . , L, with the concentration 
p being the only free parameter. This system of equations 
is easily shown to have a unique solution 



(28) 



^" l + a + ... + a^' 
where the parameter a can be determined using 



Lp. 



(29) 



The concentration p is a monotonic function of a, increas- 
ing from for a = to 1 in the limit a ^ oo. The value 
a = 1 corresponds lo p — \ and, generally, 

p{a)^\-p{\la). (30) 

Having obtained Q„ we can calculate the current as 



Using 



(ELo ^ 
it is easy to see that 

J(p) = J(l-p). 



(31) 



(32) 



Because for p ^ i Eq. ( p9[ ) implies a « pjL, using our 
formula (|3^) we conclude that for small concentrations 
of particles the current J grows linearly with p, 

J K.h-^-i^L-^p. 



(33) 



For p = \ the mean-field theory (P^) predicts 



J(0.5) 



(34) 



C. Numerical simulations 

In our simulations we used a lattice with iV = 400 sites 
consisting oi p = 100 segments, each of length L = 4. 
We used a sawtooth potential with 7 = exp(— 2) « 0.135. 
The number of particles in the system varied from M = 1 
to M — 399. We carried out our simulations for t = 
10^ Monte Carlo time steps per particle and the results 
were averaged over 10 different realizations of the process, 
which enabled us to estimate the statistical errors. 

We first present simulation results for the current at 
a fixed concentration p = 0.5, or for M = 200, as a 
function of the bias parameter b for bias to the right, 
and b^^ for bias to the left, respectively. Figure || shows 



the current J observed in simulations (symbols) together 
with a simple approximation obtained by multiplying a 
single-particle current Ji [ pT| by the number M of par- 
ticles in the system (free particle approximation). One 
observes that the current in the case of a system with a 
hard-core interactions is reduced as compared to the case 
of non- interacting particles; but the general behavior as 
a function of the bias parameter is practically the same. 
In particular, the rectification effects for particle motion 
in nonsymmetric potentials are qualitatively the same in 
both cases. The inset in Fig. ^ depicts the ratio J/MJi 
as a function of the bias. Owing to (^ij) we expect that 
for 6 > 1 J/MJi = {N - M)/{N - I) ^ I - p. For 
6 = 20 we found J/MJi w 0.5014 ± 0.0001, in excellent 
agreement with the theoretical value 200/399 w 050125. 
For 6 <IC 1 our mean-field approximation ( |3^ ) predicts 
J/MJi = LV2(L + 1)2 = 0.32; for 6 = 1/20 our simula- 
tions yielded a slightly smaller value 0.305 ± 0.001. 

We now discuss the dependence of the current on con- 
centration for selected values of the bias 6 > 1, or 6 < 1, 
respectively, and compare the results with the theoretical 
considerations of Sec. IIIB. In Fig. ^ we present results 
of our simulations for a bias to the right (6 = 30, 10 and 
2). For a strong bias {b = 30) the agreement with the 
theoretical prediction, Eq. (p^), is very good. 

The results obtained for a bias to the left (6 — 0.001, 
0.1, 0.5, and 0.9) are depicted in Fig. ||. We can see that 
if the bias is strong (6 < 0.1), the agreement between the 
mean-field theory (solid line) and the simulation data 
(circles and crosses) is very good for concentrations close 
to and 1. However, for p w i we observe that the 
mean-field theory tends to overestimate the actual value 
of J by approximately 5%, which is much more than the 
statistical errors of our data (the relative standard de- 
viation at p = 0.5 is about 0.33%). We repeated our 
simulations for larger number of Monte-Carlo time steps 
{t = 5 ■ 10^) and for different values of the bias b, but the 
difference between simulations and the theory remained 
practically the same. We thus conclude that it is not a 
numerical artifact. A similar discrepancy was observed 
by Tripathy and Barma |l^], who considered a TASEP 
with random transition rates. However, in their model 
the mean-field approach underestimated the magnitude 
of the current obtained in simulations for p « 0.5. More- 
over, they found that J(p) has quite a broad plateau 
around p = 0.5. This phenomenon is not observed in our 
case because the transition rates in our model are not 
random. 



IV. SYMMETRY PROPERTIES 

In this section we discuss the symmetry properties of 
our lattice-gas model with nonsymmetric potentials, and 
of related models. In the simulations, as well as in the 
mean-field approximation, the current exhibits a particle- 
vacancy symmetry. 
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J{p) = J(l - p) . (35) 

The symmetry properties of the TASEP have been an- 
alyzed in [p^-pO[ and the relation Eq.(^5|) has been es- 
tablished in this context. However, the model employed 
in those references differs in important aspects from our 
model. Hence a detailed discussion is in order. 

The particle-vacancy symmetry of the current for the 
TASEP has been shown in Refs. ||T^-pO|| for disordered 
hopping potentials where the transition rates are associ- 
ated with the bonds between the sites. If the motion of 
the particles is reversed (symmetry operation T accord- 
ing to Refs. jl^j2^, the particles experience the same set 
of transition rates as before, only the order of the rates 
has been changed. If the vacancies are interpreted as par- 
ticles (symmetry operation C), they experience the same 
transition rates as the particles after the operation T. 
The symmetry under CT is evident; the nontrivial state- 
ment is the symmetry of the current (up to a sign) under 
the operations C, or T, separately. 

The class of models for the hopping potential that are 
considered here do not correspond to bond disorder. The 
set of "right" transition rates is different from the set of 
the "left" transition rates. If a strong bias h » 1 to 
the right is applied, leading approximately to a TASEP, 
the current is different from the case of strong bias to 
the left with h^^ » 1. In other words, the symmetry 
under reversal of motion T does not exist for the class 
of models leading to rectification, by their definition. If 
the vacancies are considered as particles, they experience 
the same set of transition rates as the original particles, 
see also below. We conjecture that symmetry under the 
operation C also exists for our models, if the limiting case 
of the TASEP is considered. Hence we expect Eq.(|35|) 
to be approximately valid for the models that lead to 
rectification effects, in the limit of very strong bias. 

The sawtooth potential that is investigated in this pa- 
per has a special symmetry which will be described now. 
In the limit of concentration of the lattice gas approach- 
ing one, the particle problem is equivalent to the problem 
of hopping motion of single, independent vacancies. The 
hopping transitions of an isolated vacancy are reversed in 
comparison to the transitions of the particle that makes 
an exchange with the vacancy, e.g., 

^li+i = Tu+i ■ (36) 

Using the rates Eq.(^6|) it is easy to reconstruct the 
hopping potential for single vacancies. If this construc- 
tion is done for the the extended sawtooth potential of 
Fig. 1(b), a sawtooth potential is obtained for the vacancy 
which is mirror-symmetric with respect to the original 
sawtooth potential, see Fig. ^ If a bias is applied to the 
particles, expressed by the factor b in the transition rates 
to the right, the factor b appears in the transition rates of 
the vacancy to the left. It is evident from this considera- 
tion that the particle current for p — > is identical to the 
one for p — > 1. It is obvious that a particle- vacancy sym- 
metry pertains for the problem of motion of lattice gases 



in a sawtooth potential with the above symmetry prop- 
erty; hence we expect Eq.(|3^) to be valid for all values of 
the bias b. 

We point out that the sawtooth potential represents a 
special case; general nonsymmetric potentials do not pro- 
vide mirror-symmetric potentials for the vacancies in the 
limit p ^ \. For instance, if the potential corresponding 
to an Ehrlich-Schwoebel barrier (see, eg. |^) is trans- 
formed by using Eq.(|36|) into the corresponding hopping 
potential of a single vacancy, a different potential is ob- 
tained. As a consequence, the mobility of a single particle 
is different from the mobility of a single vacancy. Hence 
for this example J{p) ^ J(l — p) for b close to 1. This 
example is sufficient to show that the particle-vacancy 
symmetry (|35| ) cannot be generally valid, for arbitrary b. 
Another counterexample is provided by the random-trap 
model, see Ref.( ||). 

V. CONCLUDING REMARKS 

In this paper we investigated the motion of lattice-gas 
particles in hopping potentials that are composed of seg- 
ments without mirror-reflection symmetry. We consid- 
ered in particular the effects of exclusion of multiple oc- 
cupancy of sites, under various bias conditions. We first 
studied the case of two particles on a ring of 4 sites with 
a sawtooth potential. The explicit solution of this simple 
system can be given, and interesting conclusions emerge 
in the limits of large bias to the right, or to the left. We 
point out that the ring with 4 sites is a model case for 
the treatment of 2 site-exclusion particles on a finite ring; 
larger systems can be solved in a similar manner, e.g., by 
using symbolic formula manipulation programs. 

We then investigated the case of many particles on ex- 
tended systems which consist of periodic repetitions of 
sawtooth potentials. These systems behave, for strong 
bias in one direction, as uniform systems where the result 
for the current of lattice gases is known. For strong bias 
in the reverse directions, the extended sawtooth potential 
acts as a periodic arrangement of weak links. A mean- 
field expression for the current can be derived for this 
case from the cluster dynamics of the particles on the seg- 
ments, which shows similarities to the cluster dynamics 
of the bosonic lattice gases of Ref. . Good agreement 
with the numerical simulations was found for both cases 
under strong bias; deviations appear for smaller bias val- 
ues. The results for the current exhibit a particle- vacancy 
symmetry as a consequence of a special particle-hole sym- 
metry of the hopping processes in the sawtooth potential 
used. 

Generally, the current per particle of a site-exclusion 
lattice gas shows the same qualitative behavior as a func- 
tion of the strength and the direction of the bias pa- 
rameter, as the current of independent particles. This 
observation is important for possible applications, for in- 
stance for transport through channels in membranes or 
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through layered structures with suitable potential struc- [23] K.W. Kehr, O. Paetzold, and T. Wichmann, Phys.Lett.A 
tures. It means that qualitative or even semiquantitative 182, 135 (1993). 

predictions of the effects of strong bias on the current can [24] R. Kutner, K.W. Kehr, W. Renz, and R. Przenioslo, 
already be obtained from the single-particle description. J.Phys.A: Math. Gen. 28, 923 (1995). 
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FIG. 1. (a)Ring of 4 sites with 2 particles, (b) Sawtooth 
potential with period 4 without bias (6 — 1) with 2 transition 
rates indicated. 
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b, 1/b 

FIG. 2. Particle current J (arbitrary units) as a function 
of the dimensionless bias parameter. Upper curves: bias to 
the right, ordinate indicates b; lower curves: bias to the left, 
ordinate indicates 1/b. Dashed lines: 2 particles on the ring 
with 4 sites; Full lines: single particle on the ring. Inset: 
behavior for small bias (linear axes). 



7 




10" 



10" 



10' 



b, 1/b 



FIG. 3. Current J (arbitrary units) as a function of the 
dimensionless bias parameter for the concentration p — 0.5. 
Upper curves: bias to the right, ordinate indicates b; lower 
curves: bias to the left, ordinate indicates 1/b. Full lines: sin- 
gle-particle current Ji of Ref. multiplied by the number 
of particles M. Symbols: result of numerical simulations for 
N = 400, L = 4, 7 = exp(-2), i = 10^. The (semilogarith- 
mic) inset shows the ratio J/MJ\. 
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FIG. 5. The dimensionless current, J/&7~ , as a function 
of the dimensionless concentration p for the bias 6 = 10"'^ 
(crosses), 0.1 (circles), 0.5 (squares) and 0.9 (diamonds). The 
solid line was computed using Eq.(^l|). The parameters are 
the same as in Fig.W. The error bars, shown only for b = lO"'^, 
are of similar order for other values of b and represent the 
standard deviation. 
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FIG. 4. The dimensionless current, J/jby, as a func- 
tion of the dimensionless concentration p for the bias 6 = 30 
(crosses), 10 (circles), and 2 (squares). The solid line was com- 
puted using Eq. (p2|) . The parameters are A'' = 400, L = 4, 
7 = exp(— 2), t = 10®. Results were averaged over 10 MC 
simulations. The error bars are shown only for b = 30; for 
other values of b they are similarly small. 
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FIG. 6. (a) Repetition of the sawtooth potential of Fig.l 
with lattice-gas particles indicated, (b) Effective potential for 
a single vacancy, constructed according to Eq.Bg. 
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